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ABSTRACT 

We present three-dimensional numerical simulations of a magnetic loop evolving in 
either a convectively stable or unstable rotating shell. The magnetic loop is introduced 
in the shell in such a way that it is buoyant only in a certain portion in longitude, thus 
creating an fi-loop. Due to the action of magnetic buoyancy, the loop rises and develops 
asymmetries between its leading and following legs, creating emerging bipolar regions 
whose characteristics are similar to the ones of observed spots at the solar surface. 
In particular, we self-consistently reproduce the creation of tongues around the spot 
polarities, which can be strongly affected by convection. We moreover emphasize the 
presence of ring-shaped magnetic structures around our simulated emerging regions, 
which we call "magnetic necklace" and which were seen in a number of observations 
without being reported as of today. We show that those necklaces are markers of 
vorticity generation at the periphery and below the rising magnetic loop. We also find 
that the asymmetry between the two legs of the loop is crucially dependent on the initial 
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magnetic field strength. The tilt angle of the emerging regions is also studied in the 
stable and unstable cases and seems to be affected both by the convective motions and 
the presence of a differential rotation in the convective cases. 

Subject headings: convection, magnetic fields, magnetohydrodynamics (MHD), meth- 
ods: numerical. Sun: interior. Sun: rotation 



1. INTRODUCTION 



Our turbulent active Sun possesses a rather regular magnetic cycle that has been observed di- 
rectly through the photospheric emergence of sunspots for about four centuries now. This magnetic 
flux emergence, which occurs daily at the Sun's surface, appears not only at large scales (sunspots 
and active regions) but also at scales as small as the intra-network field and there is indication that 
even smaller scale flux emerges at the solar surface, which remains unresolved with the present 
magnetographs (Sanchez Almeida 1998; van Driel-GeszteIyil|2002D . However, since a huge amount 
of flux is brought toward the photosphere by solar active regions and because of the well-established 
robust properties of such large-scale emergence, we focus our studies on the understanding of the 
dynamical evolution of large and strong magnetic structures rising from the solar interior to the 
surface. The appearance of active regions can be summarized in a few steps: small magnetic bipoles 
first emerge all over the future active region and then merge to form sunspots at the outer edges 
of the emerging flux region (EFR). The growing phase of a typical EFR lasts about 3-5 days. The 
flux in an active region is usually larger than lO^^Mx in each polarity. The orientation of the bipole 
may be arbitrary at first, but generally in 1-3 days it becomes slightly tilted (3^-10^) with respect 
to the east-west direction, with the leading polarity closer to the equator, thus conforming to Joy's 



law (Hale et al. 1919). The tilt angles of individual active regions show a large scatter about the 
mean but the measurements indicate a tendency for the mean value to linearly increase with the 
latitude of emergence dWang fc SheeIeylp89l[l99T||D'Silva fc Choudhuri|[l993 l) . Moreover, various 
asymmetries have been detected in EFRs in particular in the morphology of the leading and trailing 
spots. Indeed, the leading polarity generally appears more intense and concentrated than the trail- 
ing one which tends to be more fragmented. Another interesting feature of EFRs is their inherent 



twist. Indeed, careful multi- wavelength analysis of flux emergence in active regions (e.g. Leka et al 



1996) indicates that they appear with a certain amount of twist (i.e., a non-zero current along their 
This twist is then thought to be responsible for the characteristic pattern known as tongues 



axis 



observed in newly emerging regions ( Lopez Fuentes et al.||2000 ) and identified as a major ingredient 
in mass ejections or flares at the solar photosphere (e.g.. Fan fc Gibson||200"4 ; Torok fc Kliem||2005 



Aulanier et al.|20T0 ). Other structures such as an annular shape of the magnetic field around active 
regions have been observed through various instruments (e.g., Liu fc Zhang||2006] ) but never really 
reported as identified structures. However, as we will see, they can serve as strong constraints on 
models and simulations and the relevant physical processes responsible for their existence will be 
pointed out. 
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With all those observational constraints in mind, theoreticians and modelers have tried to 
analyze in detail the dynamical evolution of strong magnetic structures from the base of the con- 
vection zone (CZ) where they are created to the surface where they appear as sunspots. A possible 
scenario for such an evolution would be that strong toroidal fields are created at the base of the 
CZ by the differential rotation shearing a preexisting poloidal field in the tachocline (see landmark 
papers of |Moffatt|1978l |Parker[T993l recent simulations of[B rownms: et al.|2006| and a review on the 
solar dynamo theory by Ossendrijver 2003 for example). After a certain period of storage in the 
subadiabatic layer below the CZ (iMoreno-Insertis et al.||1992l iSchiissler et al.||1994|), this toroidal 



field becomes subject to magnetic buoyancy instabilities (Parker 1955; Acheson 1979; Cattaneo &; 



Hughes 1988; Matthews et al. 1995) and organizes itself in buoyant arched magnetic flux tubes 



(or 0-loops, Zwann 1987). The rise of those loops under the influence of Coriohs forces, turbu- 
lent convection and large-scale flows (rotation and meridional circulation) would then produce the 
wide variety of large-scale flux emergence observed at the solar surface, with the robust properties 
summarized above. Several numerical experiments of this last step have been performed using 
various levels of approximation, from the classical "thin flux tube" models ( Spruit||1981 ; Choudhuri 



& Gilman 1987) to the more recent two-dimensional and three-dimensional calculations (see Fan 



2004 for a detailed review on this subject). In particular, Emonet &; Moreno-Insertis (1998) showed 



the necessity to start with a twisted magnetic field to counteract the vorticity generation within 
the flux tube and its subsequent breakup. Fan et al. (2003) demonstrated the significant effect of 
convective motions on the deformation of the magnetic loop when the initial field was comparable 
to the equipartition field strength (corresponding to a magnetic energy equals to the kinetic energy 
of the strongest downflows). More recently, Jouve &; Brun (2009) considered the effects of spherical 
geometry, convective motions, and self-consistently developed mean flows (differential rotation and 
meridional flow) on a rising flux tube. It was found that an initial twist of the field lines of about 
14 turns along the 360^ longitude as well as an initial field strength of about 130kG were necessary 
for a coherent radial rise of the tube, that the tilt angle was significantly influenced not only by the 
initial twist and Coriolis force but also by the local convection, and that meridional circulation was 
able to advect initially weak flux tubes toward the poles when they reached the top of the domain. 
They showed moreover a quantitative study of the effect of magnetic diffusivity by varying the 
Prandtl number in their simulations. Fan (2008) also investigated the effects of spherical geometry 
but on an fi-loop (a flux tube made buoyant in a limited range of longitudes) and thus focused on 
the properties of simulated individual active regions, but without convection. We here propose to 
compute the dynamical evolution of an fi-loop both in an isentropic and fully convective spherical 
shell, seeking to answer the following questions. 



- What are the physical processes at the origin of the observed tilt angles of active regions and 
is it possible to reproduce Joy's law? 

- What is the influence of convective motions on a rising f]-loop and the subsequent bipolar 
emerging regions? 

- Can we quantify the effects of large-scale mean flows? 
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How important is the twist of the magnetic structures in the evolution within the convection 
zone and on the morphology of EFRs? 

How does the structure of simulated active regions compare to observations? 



This work is thus an attempt to answer those questions using numerical simulations of a 
convective spherical shell, using the Anelastic Spherical Harmonic (ASH) code which solves the 
MHD equations under the anelastic approximation, valid within the bulk of the convection zone. 
We do not address here the evolution of emerging flux in the solar atmosphere which needs to be 
computed in a fully compressible simulation. Extensive studies of magnetic flux emergence in the 



atmosphere have been performed (e.g. Magara 2004; Manchester et al. 2004; Murray et al. 2006 



Martinez-Sykora et al. 2008; Archontis &; Torok 2008) which showed that magnetic flux reaching 



the photosphere could undergo a dynamical expansion into the atmosphere as a result of the non- 
linear growth of the magnetic buoyancy instability and that it is difficult for a twisted flux tube 
to rise bodily into a corona because of the heavy plasma trapped in the magnetic dips. Successive 
magnetic reconnections could however help to achieve actual emergence of twisted magnetic loops 
into the corona dPariat et aL][2004 |. 

This article is organized as follows. Section [2] presents the model chosen, the hydrodynamic 
background as well as the initial conditions which will favor the emergence of individual active 
regions. The following short Section [3] presents the various cases calculated, which will be consid- 
ered in an isentropic background (without convection) and with convective motions acting on the 
magnetic loop. Sections |4] and [5] respectively presents the results of the MHD simulations in the 
isentropic and convective cases. Finally, we discuss the results and conclude in Section [6j 



2. THE MODEL 

2.1. Anelastic MHD Equations 

All calculations presented in this work were computed with the ASH code, which solves 
anelastic magnetohydrodynamic equations in a rotating spherical shell. ASH uses a pseudo-spectral 
method in space and semi-implicit approach in time (e.g., Clune et al.|1999 ; Miesch et al.|2000 ; Brun 



et al.|[2004 ). The full nonlinear evolution of the velocity and magnetic fields is computed, whereas 
the equations are linearized in thermodynamic variables with respect to a spherically symmetric 
mean state to have density p, pressure P, temperature f , and specific entropy S. Perturbations 
are denoted as p, P, T, and S. The equations being solved are 



V • (pv) = 0, 
V-B = 0, 



(1) 
(2) 
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dv 

+ (v • V)v + 2!^o X v] = - VP + pg 



,dS 



(3) 



+ — (VxB) xB-V-P-[VP-pg], 
47r 



pT— + pTv • V(5 + 5) = V • [KrpcpV{T + T) 



(4) 



+ KopfVS + KpfVS] + ^f 



+ 2pv 
d'B 



dt 



= Vx(v X B) - Vx(r7V x B) 



(5) 



where v = (vr^ve^Vcj)) is the local velocity in spherical coordinates in the frame rotating at a 
constant angular velocity Qq, g is the gravitational acceleration, B = (B^, Bq, B^f^) is the magnetic 
field, j = (c/47r)(V x B) is the current density, Cp is the specific heat at constant pressure, /^^ is 



the radiative diffusivity (fitted to a one-dimensional seismically calibrated solar model, see Brun 



77 is the effective magnetic diffusivity, and V is the viscous stress tensor. To represent 
the unresolved subgrid-scale processes, the code uses an effective eddy viscosity u and an effective 
thermal eddy diffusivity. They are both here to model the enhanced diffusion associated with 
turbulent convective motions in the bulk of the convection zone. In the present simulations, they are 
chosen to vary in radius only as p~^l^ . The thermal diffusion /^o acting on the mean entropy gradient 
(contrary to the thermal eddy diffusivity n which acts on the entropy perturbation) occupies a 
narrow region in the upper convection zone. Its purpose is to transport heat through the outer 



surface where radial convective motions vanish (Oilman &; Glatzmaier 1981; Wong &; Lilly 1994) 



and is thus chosen such that the entire luminosity will be transported at the top of the domain by 
the so-called unresolved heat flux Fu — —K^pTdSjdT. The diffusivity k is purely dissipative and 
acts to smooth out entropy variations, whereas a^o is essentially a cooling term near the top of our 
computational domain. To complete the set of equations, we use the linearized equation of state. 



P_ 
P 



P 
P 



T 
T 



P 



(6) 



where 7 is the adiabatic exponent (equal to 5/3 here), and assume the ideal gas law 



p = npf, (7) 

where TZ is the ideal gas constant, taking into account the mean molecular weight p corresponding 
to a mixture composed roughly of 3/4 of hydrogen and 1/4 of helium per mass. The reference or 



mean state (indicated by overbars) is derived from a one-dimensional solar structure model (Brun 



et al. 2002). It begins in hydrostatic balance so that the bracketed term on the right-hand side of 
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Equatioi][4] initially vanishes. However, as the simulation evolves, turbulent and magnetic pressures 
drive the reference state slightly away from strict hydrostatic balance. 

The ASH code uses a spectral decomposition for all the variables, namely, spherical harmonics 
in the horizontal direction and Chebyshev polynomials in radius. The collocation points are thus 
the zeros of the Legendre polynomials in latitude, uniformly distributed in longitude, and the 
distribution is denser at both ends of the domain in the radial direction. The code integrates the 
MHD equations presented above over O.72i?0 < r < 0.96Rq (with an overall density contrast of 
about 24 between top and bottom), < < tv and < < 27v. The resolution is chosen here such 
that at least 20 points in radius and in latitude are contained inside the flux tube section. The 
minimal resolution used here is thus Nr x Nq x N^fy = 256 x 512 x 1024. 

Finally, the boundary conditions for the velocity are impenetrable and stress-free at the top 
and bottom of the shell. We impose a constant entropy gradient top and bottom for the isentropic 
case and for the fully convective case, a latitudinal entropy gradient is imposed at the bottom, as 



in.Miesch et al. (2006). In all cases, we match the magnetic field to an external potential magnetic 



field at the top and the bottom of the shell (Brun et al.||2004). 



2.2. Initial Conditions to Favor the Creation of fi-loops 

In this paper, we implement initial conditions for the magnetic field similar to what was 



used in Jouve &; Brun (2009). An initially axisymmetric magnetic structure is embedded in an 
unmagnetized stratified medium. In order to keep a divergenceless magnetic field, we use a toroidal- 
poloidal decomposition. 



B = V X Vx(Ce^) + Vx{Aer) 
the expressions used for the potentials A and C for the flux tubes are: 



(8) 



A = —Aq r exp 



C = -Ao^q exp 



r-Rt 



r-Rt 



1 + tanh 2 



a/Rt 



1 + tanh 2 



a/Rt 



(9) 



(10) 



where Aq is a measure of the initial field strength, a is the tube radius, (Rt^Ot) is the position of 
the tube center and q is the twist parameter. 



In Jouve & Brun (2009), we derived an expression for the winding degree of the field lines (i.e., 
the number of turns that the field lines make over the whole tube length 27vRtsm0t, assuming a 
uniform twist): 
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TvRt sm 9t , 

n = tan ip, 

2a 



(11) 



where the tangent of the pitch angle t/j is related to the twist parameter via the following 
expression, considering that we are 6 — 6t and at the tube periphery r = i?^ + a, as in | Jouve"fc 



Brun (2009): 



tan^ 



qa 



Rt + a 



(12) 



In all cases, the tube radius is set to a = 2 x 10^ cm, about a 10th of the depth of the modeled 
convection zone and is introduced at the base of the CZ at Rt = 5.2 x 10^^ cm. If we consider that 
the typical flux in an observed active region is 10^^ — lO^^Mx and that an initial flux tube has a 
magnetic field of about 5 x lO^G as in our simulations, the conservation of flux tells us that the 
radius of the loop should be in reality closer to 8 x lO^cm. Our flux tubes are thus bigger than 
what could be expected. Different values for the loop radius have been shown to influence flux 
emergence ( D'Silva fc Choudhuri||1993 ; Emonet &; Moreno-Insertis 1998, see also the discussion of 
Section [g]). We here limit our parameter study to a fixed loop radius (see Jouve & Brun 2009, for 
a different set of parameters). The initial field strength Aq, the initial twist of the field lines q as 
well as the latitude of introduction Of will be varied in our models to investigate the influence of 
these various parameters. 

In order to get a flux tube buoyant on a small portion in longitude only, we initially apply a 
perturbation on the background entropy field located at the initial position of he flux tube in (r, 0) 
and possessing a Gaussian profile in longitude 0. The entropy perturbation has thus the following 
expression: 



Sin S'ext 



As X exp 
2Rt 



r-Rt 



cosh^ 



/Rt) 



(13) 



C^ + exp -(^^— ^) 



where As is the amplitude of the entropy perturbation, Cs is a constant value controlling the 
buoyancy of the rest of the loop, 00 is the longitude of maximum perturbation (or buoyancy), and 
(pe is the extent of the entropy perturbation. 



The effect of such a perturbation on entropy is to produce an additional density deficit inside the 
flux tube at a particular location in longitude. As a consequence, we can derive the new maximum 
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density contrast between the tube and its surroundings, which wih result from this additional 
perturbation. Total pressure equilibrium is not enforced at the beginning of the simulation and 
since an initial magnetic tension is also present in most of our cases (the twisted cases), the loop 
is not initially in mechanical equilibrium outside the portion which is made buoyant, even if Cs is 
modified. However, there is an initial adjustment visible at the very early stages of the dynamical 
evolution where the pressure perturbation inside the loop will tend to compensate for the magnetic 
pressure, in particular the loop gets slightly squeezed in an Alfven time at the very beginning of 
the simulation, without much influence on the following evolution. We thus derive here only an 
indication for the efliciency of the magnetic buoyancy in a situation of pressure equilibrium and 
entropy perturbation (not exactly satisfled here) in writing the following relations respectively for 
the total pressure and the entropy: 



PL.-PL-^ (14) 



Max (Sin - Sext) = AS = As X X (1 - Cs) , (15) 

a 

where P^^^ and P^^ are the gas pressure respectively outside and inside the tube and where the 
maximum of iSin — S^ext is given by Equation [13] and taken r — Rt^ 9 = 6t^ and = ^q. We then 
use an additional relation to make the density p appear: 



S — CyhiP — Cp\n (16) 

with Cp the speciflc heat at constant pressure (equal here to 3.4 x 10^erg.K~^.g~^), the speciflc 
heat at constant volume, and 7 = Cp/c^ = 5/3 the adiabatic index. 

These equalities lead to the following relation between entropy, pressure, and density inside 
and outside the flux tube: 

A5 = Max (Sin - Sext) = Max (c^ In ^ - Cp In ^\ (17) 
and thus we flnd an expression for the maximum buoyancy of our magnetic structures: 



P \ Pext J V Cp y V SvrPext / \ '^P J \ 87r7Pext / 

(18) 

where Bq will be the maximum fleld strength on the tube axis. A flrst-order Taylor expansion 
has been done to give an approximate value of the buoyancy, considering that both AS/cp and 
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BQ/87rPext are equally small values. In the remaining of the paper, we will use the following 
parameters for the standard cases for the entropy perturbation: Cs — 0.35, (j)e — 15^, 00 = 100^ and 
the average external pressure used to estimate the maximum buoyancy will be Pgxt = 4 x lO^^dynes. 
We use a fixed value for Cs since we do not seek to maintain a perfect mechanical equilibrium for 
the rest of the loop, otherwise the value of Cs should depend on the magnetic field. 



2.3. The Background Hydrodynamical Models 



The same hydrodynamical background was used in this work as in Jouve &; Brun (2009). 
Namely, two situations will be investigated. In the first, the magnetic structure is introduced in 
an isentropic layer with solid body rotation and without convection. In the second situation, the 
convective instability is triggered and differential rotation as well as a large-scale meridional flow are 
self-consistently created within the bulk of the convection zone. The competition between buoyant 
loops and convective motions and large-scale flows will in this case be studied and compared to 
the stably stratified situation. In particular, we focus on the shape of the developing loop, its 
asymmetries and the morphology of the emerging radial field it generates. Moreover, we focus 
in this study on the comparison to observations of large-scale EFRs, extrapolating the magnetic 
configurations we get from our simulations at the upper boundary (i.e., O.OGi?©) to what could 
happen higher up, closer to the photosphere. We specifically investigate the evolution of flux, tilt 
angle and spot morphology during emergence. 



3. CASES CALCULATED 

We investigate in this study seven different cases. All parameters used for each case are 
summarized in Table [ij Four of the cases go by pair with the same buoyancy (but different field 
strength and entropy perturbation: Cases 1 & 3 and 2 &; 4) and the last three are extreme cases (the 
weakest Cases Oa and Ob and strongest field Case 5). For Cases 1, 2, 3, and 4, the twist parameter 
was also varied. The cases with a right-handed (or Positive) twist (corresponding io q — 30) will be 
called by the case number with the suffix TwP, and for the left-handed (or Negative) twist, we will 
add TwN. For example, the case where = 5 x lO^G, AS = 3380 (corresponding to a buoyancy 



Ap/p of 1.1 X 10 ^ according to Equation 18) and which has a right-handed twist g = 30 will be 
called Case ITwP. 

Some additional cases were also computed: 

Case 1 with the loop introduced at the latitude of 60^: Case lTwP60 and 15^: Case lTwP15. 
Case 3 without twist: Case 3NoTw and with the twist at its threshold value (defined in the 



last paragraph of Section 4.1.1): Case 3TwTh 



The values of the entropy perturbation were chosen to ensure a rise time of the same order 
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Table 1. Key Parameters of the Various Cases. The amphtude of the magnetic field is given in 





latitude in 


degrees and entropy perturbation 


iXo m ' 


erg.K 












Bo 


Sg/87rPext 


Twist q 


Latitude 


AS 


AS/cp 


Buoyancy Ap/p 


CaseOa 


2 X 10"^ 


4 X 10" 


7 


30 


30 


1690 


5 X 


10-6 


5.1 X 


10" 


-6 


CaseOb 


4 X 10"^ 


1.6 X 10" 


-6 


30 


30 


1690 


5 X 


10-6 


5.9 X 


10" 


-6 


Case ITwP 


5 X 10^ 


2.5 X 10" 


-6 


30 


30 


3380 


9.9 > 


< 10-6 


1.1 X 


10" 


-5 


Case ITwPlS 


5 X 10^ 


2.5 X 10" 


-6 


30 


15 


3380 


9.9 > 


< 10-6 


10 


-5 




Case lTwP60 


5 X 10^ 


2.5 X 10" 


-6 


30 


60 


3380 


9.9 X 10-6 


10 


-5 




Case ITwN 


5 X 10^ 


2.5 X 10" 


-6 


-30 


30 


3380 


9.9 > 


< 10-6 


10 


-5 




Case 2TwP 


10^ 


9.9 X 10" 


-6 


30 


30 


3380 


9.9 > 


< 10-6 


1.6 X 


10" 


-5 


Case 2TwP15 


10^ 


9.9 X 10" 


-6 


30 


15 


3380 


9.9 > 


< 10-6 


1.6 X 


10" 


-5 


Case 2TwP60 


10^ 


9.9 X 10" 


-6 


30 


60 


3380 


9.9 X 10-^ 


1.6 X 


10" 


-5 


Case 2TwN 


10^ 


9.9 X 10" 


-6 


-30 


30 


3380 


9.9 X 10-6 


1.6 X 


10" 


-5 


Case 3TwP 


10^ 


9.9 X 10" 


-6 


30 


30 


1690 


5 X 


10-6 


10 


-5 




Case 3NoTw 


10^ 


9.9 X 10" 


-6 





30 


1690 


5 X 


10-6 


1.1 X 


10" 


-5 


Case 3TwTh 


10^ 


9.9 X 10" 


-6 


10 


30 


1690 


5 X 


10-6 


1.1 X 


10" 


-5 


Case 3TwN 


10^ 


9.9 X 10" 


-6 


-30 


30 


1690 


5 X 


10-6 


1.1 X 


10" 


-5 


Case 4TwP 


1.35 X 10^ 


1.8 X 10" 


-5 


30 


30 


1690 


5 X 


10-6 


1.6 X 


10" 


-5 


Case 5TwP 


2 X 10^ 


4. X 10- 


5 


30 


30 


3380 


9.9 > 


< 10-6 


3.4 X 


10" 


-5 
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of the convective turnover time, i.e., between 10 and 40 days approximately. We can compare the 



buoyancy of the loops with and without the entropy perturbation by setting AS = in Equation 18 



We find that the value of Ap/p is multiplied by 7.65 when the entropy perturbation is introduced 
in Case 1 (compared to the case without entropy perturbation) and by a factor 1.41 in Case 5, 
so that the stronger field cases are less influenced by the entropy perturbation applied. Moreover, 
the values of the entropy perturbations used here are comparable in amplitude to the typical 
fluctuations measured at the base of our computational domain in our convective runs. 



4. EVOLUTION OF A SINGLE LOOP WITHOUT CONVECTION 

We first compute the evolution of a single fi-loop evolving in a stable background with respect 
to convection. It will first confirm that we indeed correctly control the buoyancy of the flux tube 
with the initial conditions we apply. Secondly, it will help us see if the important parameters for 
the evolution of a uniformly buoyant flux tube are still the same here and what differences may be 
pointed out. 



4.1. Favorable Conditions for Emergence 

4.1.1. Twist 

We know from previous experiments that a certain amount of twist is needed for a magnetic 
flux tube to be able to rise coherently in a convectively stable (or unstable) background (e.g.. 



Emonet &; Moreno-Insertis 1998). Indeed, in a non-twisted case, the azimuthal vorticity created 
by the gravitational torque acting on the loop cannot be compensated. Two counter vortices 
thus quickly appear and break the loop apart, preventing it from emerging at the surface. It 
was argued previously that this situation may differ when an fi-loop is considered, instead of a 



uniformly buoyant structure (Abbett et al. 2000). The additional magnetic tension coming from 



the deformation of the magnetic torus may oppose the vorticity generation. 

Figure [l] shows the section of an untwisted fi-loop corresponding to Case 3NoTw of Table 
[Tj The zoom is made on the longitude where the loop is caused to be buoyant, thanks to the 
combined effects of magnetic pressure and entropy perturbation. Two instants in the evolution 
are shown. We represent here the contours of the azimuthal (or axial) component of the magnetic 
field superimposed on the azimuthal vorticity. It is rather clear in this case (especially on the left 
panel) that strong vortices are formed inside the flux tube, which have the tendency to significantly 
deform the magnetic structure. However, as the loop rises further up, it becomes more arched and 
more magnetic tension is generated. As a consequence, even without any twist, some coherence is 
maintained in this case, even if a large amount of flux is lost in the process. 

This flux loss is visible in Figure [2] where the profile of the toroidal magnetic field is shown at 
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Time=5 days Time=13.6 days 



Fig. 1. — Cross section of the untwisted fi-loop of Case 3NoTw, introduced at a latitude of 30^, 
at two instants in the evolution, 5 and 13.6 days. Magenta contours of the azimuthal component 
of the magnetic field are superimposed to the azimuthal vorticity which presents a characteristic 
profile of two counter vortices. 



4x10* 




40 50 60 70 80 

Colatitude (Degrees) 



Fig. 2. — Profile of the toroidal field at the radius where it peaks after 6 days of evolution (i.e., 
O.82i?0). The three cases represented are the untwisted case {q — 0) Case 3NoTw, the twisted 
case Case 3TwP {q — 30), and the intermediate case Case 3TwTh where the twist value is at its 
threshold {q = 10). 

the radius where it peaks after about 6 days of evolution, as a function of colatitude. We remind 
that the flux ropes were introduced here at 60^ of colatitude and with an amplitude of lO^G. In 
the non- twisted case (dashed line on the figure), we clearly see that the two counter vortices have 
created two distinct magnetic flux concentrations, whose maximum amplitude reaches a value of 
2.5 X 10"^ G while the maximum amplitude of the twisted case is 4 x lO^G. The flux rope has been 
eroded by the vorticity generation and the two vortices tend to dominate the further evolution 
of the magnetic structure which separates into two flux concentrations. However, along the loop 
(following its axis), this double-peak structure is not as clearly visible as at the longitude we chose 
here. This feature may be the consequence of the asymmetry which develops between the leading 



- 13 - 



and trailing leg. We will come back on this in the next section. 



We decided to define the twist threshold as the value above which the double-peak structure is 
not present anymore at any longitude. With this criterion, we find that the twist threshold for this 
case is g = 10, corresponding to 0.54 turns along the 15^ of extension of the fi-loop. We note that 
this threshold value is reduced by a factor 1.8 compared to the uniformly buoyant flux tubes of 



Jouve &; Brun (2009). This is due to the additional tension coming from the arching of the fi-loop 



in this case, which has already the effect of limiting the vorticity generation inside the magnetic 
structure. This result is in agreement with the work of Abbett et al. ( 2000| ). The case where the 
twist is set at its threshold value is also represented in Figure [2] (dotted line). In this case, the 
maximum value of the field is reduced compared to the more twisted case but the coherence of the 
loop is kept. 



4.1-2. Buoyancy versus Rotation 



The buoyancy force is proportional to Ap which is given by expression (18). While the loop 
rises because of the buoyancy force, a retrograde longitudinal flow develops in order to conserve 
angular momentum (see section below). As a consequence, a Coriolis force oriented toward the 
rotation axis and perpendicular to this axis also develops, opposing the component of the buoyancy 
force perpendicular to the rotation axis. The only uncompensated force is then parallel to the 
rotation axis and oriented toward the pole. If the buoyancy force is not strong enough, the loop 
will thus have the tendency to rise parallel to the rotation axis and emerge at a much higher latitude 
than where it was introduced. This is the case at fixed AS when the magnetic energy is not strong 
enough. We recover this situation when the poloidal field is too weak (corresponding to a low twist 
value, enough to maintain coherence of the loop but not enough to make it rise radially). 

Figure [3] is an illustration of a case where the rise is mainly parallel to the rotation axis. In 
this Case Oa, the twist was set to g = 30 but the initial axial field strength was set to 2 x lO^G, 
and the entropy perturbation to the rather weak value of A^- = 1680, corresponding to a density 
perturbation of Ap/p = 5 x 10~^. In this case, after 15 days of evolution, the trajectory of the loop 
is clearly parallel to the rotation axis, since the Coriolis force has strongly acted on the magnetic 



field. Indeed, as shown by Choudhuri &; Oilman (1987) in the framework of the thin flux tube 



approximation, while the loop rises, it develops a retrograde longitudinal flow in its interior. A 
Coriolis force is associated with this longitudinal flow which tends to oppose the outward component 
of the buoyancy force perpendicular to the rotation axis, leaving only the component parallel to 



the rotation axis to act on the ring. Interestingly, [Choudhuri fc Oilman (1987) find that in the 



adiabatic cases (similar to what we are studying in this section), the transition between vertical 
and radial rise happens between Ap/p = 2 x 10~^ and Ap/p = 2 x 10~^. We also find here that 
Cases Oa and Ob, where Ap/p is around 5 — 6 x 10~^, tend to be strongly influenced by the Coriolis 
force while the stronger cases, where Ap/p > 10~^, exhibit a much more radial trajectory. 
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Fig. 3. — Section of at t = 15 days at the longitude of maximum buoyancy for Case Oa. The 
trajectory of the loop is here clearly parallel to the rotation axis. 

We note that due to the low value of the buoyancy used for Case Oa, the rise of the loop is 
rather slow (about 40 days are needed to reach the top of our domain). In order to ensure an 
evolution dominated by the dynamical effects and not only by diffusion, we reduced the magnetic 
diffusivity from 7.95 x lO^^cm^.s"^ to 1.59 x 10^^, and thus used a magnetic Prandtl of 5 instead 
of 1 as with all the other cases. It does not change the fact that the rise is parallel to the rotation 
axis in this case, it is just to compute a case where the rise time is significantly lower than the 
diffusive time (which would be of order B? /rj — 58 days for the Pm — 1 case and thus too close to 
the rise time of 40 days). 

The same kind of evolution parallel to the axis of rotation is found when stronger tubes are 
introduced with a lower twist. Indeed, if the twist is low as we saw before, the amplitude of the 
magnetic field inside the loop quickly decreases until it reaches a value where the buoyancy force 
is not strong enough to keep the radial trajectory. This is visible on the right panel of Figure [l] 
where the initial axial field was lO^G but the initial poloidal field was 0. To correctly define a 
minimum value for a radial rise here, it is thus necessary to take into account the three components 
of the magnetic field rather than focus on the axial field. Let us consider a case where the entropy 
perturbation is the weakest, corresponding to a slow rise and then to a maximum effect of the 
Coriolis force. We then find that for a twist parameter g = 30 (about 1 turn over the extension 
chosen here), the initial magnetic field needed for a radial rise should be at least of 5 x lO'^G. It 
can be argued that this value is strong compared to what could be expected at the base of the 
convection zone from dynamo action but it could be reduced if a stronger entropy perturbation 
was considered, which is locally possible in the turbulent environment existing at the base of the 
convection zone. 
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To summarize the two previous subsections, we can say that to maintain a coherent radial rise 
of our fi-loop in our simulation, an initial twist of 0.54 turns along the extension of the loop, as 
well as an initial field strength of 50kG are necessary. We note that those two values are reduced 
compared to the uniformly buoyant calculations. This is promising since the observed twist values 



are of the order of less than 1 turn across the active region ( |Chae fc Moon||2005t [Demoulin et al. 



2002 ) and typical fields of more than 50kG are difficult to produce in present dynamo calculations 



dBrun et al.][2004l). 



4.2. Asymmetry between the Trailing and Leading Legs 

We now study the evolution of a sufficiently twisted and strong loop so that it will maintain 
its coherence while it rises radially. The loop is introduced at 30^ in latitude, with a maximum 
buoyancy (i.e., entropy perturbation) located at 100^ in longitude, with a twist parameter g = 30 
representing about 39 turns around the tube axis over the whole 360^ and thus corresponding to 
about 1.6 turns over the extension of the buoyant part {(j)e = 15^). With this twist parameter (as 
we saw before), the loop rises radially for a field strength of 5 x lO^G. Figure [i] shows a volume 
rendering of the toroidal field evolution in Case ITwP. 

The asymmetry between the trailing and leading legs develops during the evolution and is 
clearly seen in the last two panels of the figure, at t = 10.7 days and t = 14.7 days. This section 
will be devoted to the study of this asymmetry and the influence of buoyancy and field strength. 
The action of rotation and the associated Coriolis force over the whole loop will also be detailed 
here. This process will lead to the appearance of a tilt of the axis with respect to the east-west 
direction and will also depend on the magnetic field strength, as we will see in the following section. 

Let us enter now into the details of the asymmetry between the trailing and leading legs and the 
effect of the field strength on it. This asymmetry, already seen in thin flux tubes calculations (see 
Caligari et al.||199^ ), is due to the difference in azimuthal velocity between the top of the loop and 



its feet. Indeed, since the flux tube needs to conserve its angular momentum rsinO (rsinOQ + v^f^), 
a retrograde flow with respect to the rotating frame will develop inside the loop while it rises 
(i.e., while r increases). Since the apex of the loop is higher than the feet which stay attached to 
the bottom of the convection zone, a stronger retrograde flow appears at the apex. The profile 
of the flow along the magnetic structure then leads to the inclination asymmetry of the trailing 
and leading legs. In some of our cases, since the parameter Cs is always kept the same, parts of 
the loop located outside the buoyant portion can be slightly denser than their surroundings. As a 
consequence, those regions may sink toward the base of the convection zone and favor the difference 
in azimuthal velocity causing the asymmetry. However, in most cases (except Cases Oa, Ob, ITwP, 
and ITwN), the magnetic buoyancy can overcome this effect and the rest of the loop is slightly 
positively buoyant, thus on the contrary reducing the asymmetry. 
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Fig. 4. — Volume rendering of B^j) while the loop rises through the isentropic layer in Case ITwP. 
The loop is viewed up to the polar axis, from the South Pole. The asymmetry between the legs 
of the loop are obvious in these snapshots. The visualization was performed with SDvision@cea 
dPomarede fc Brun||20Tol ) . 

However, this explanation is valid as long as the Lorentz forces do not dominate the flow 
evolution inside the loop. Figure [5] shows a snapshot of the magnetic energy and the azimuthal 
velocity for four different cases: Cases Ob, ITwP, 2TwP, and 5TwP. We note here that Case Ob 
has a magnetic field strength of 4 x lO^G and is thus still significantly deflected toward the pole 
but is an interesting case because it maximizes the asymmetry between the legs of the loop. The 
snapshots were chosen so that all loops have reached the same height (namely O.OSi?©), which 
corresponds to different physical times since the buoyancy force (i.e., the Ap) is different in all 
cases. In order to quantify the asymmetry, we measure the angle between the tangent at the apex 
and the line going through the middle of the two footpoints and the apex. Those lines are shown in 
magenta in the top panels of Figure [5| We find that the angles have the following values: 45^, 50^, 
75^, and almost 90^ (meaning no asymmetry in this Case 5TwP where the magnetic field and the 
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Fig. 5. — Cut at the latitude of 30^ of the magnetic energy (upper panel) and azimuthal velocity 
together with the white contours of magnetic energy (lower panel) for four different cases, Cases 
Ob, ITwP, 2TwP, and 5TwP, when the fi-loops have approximately reached the same height r = 
O.93i?0. The loops are viewed down from the North Pole. The color scale used for the magnetic 
energy is the same for the first two panels and the same for the last two. 

buoyancy force are the strongest). We thus note that the asymmetry strongly decreases from Case 
ITwP to Case 2TwP, even if the density deficits are initially not very different. The reason for this 
jump in the asymmetry is the difference in field strength (lO^G for Case 2TwP against 5 x lO^G for 
Case ITwP). Indeed, in this case, the Lorentz force starts to become important in the generation 
of azimuthal velocity: strong values of Vfj) appear at the location of strong gradients of magnetic 
field. This is much more obvious in the stronger Case 5TwP. The profile of the azimuthal velocity 
shown on the last snapshot of the lower panel is clearly related to the gradients of magnetic field 
and the retrograde flow due to the conservation of angular momentum is so weak compared to the 
magnetic effects that it is not visible anymore. Instead, an azimuthal velocity directed toward the 
bisector of the tube appears on both sides of the loop (indicated by black arrows in the last panel), 
resulting in a symmetric arched structure (visible in the top right snapshot) with both feet showing 
an almost vertical orientation. 

If we write down the expression of the Lorentz force involved in the evolution of Vfj)^ we have: 



— ^ h non-magnetic terms (19j 

Be ( d dBe\ 
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This equation shows that the regions of strong magnetic field gradients wih correspond to 
the regions of production of v^^^ which can be verified from the last panel of Figure [5| where the 
concentrations of v^^ are located mainly at the footpoints (where the gradients of B^ and Bq with 
respect to are strong) and higher up at the loop periphery (where the gradients of B^^ with respect 
to r and d are strong). 

To be more precise, we can follow the evolution of the maximum and minimum of the azimuthal 
velocity which develops during the rise of the flux tube from the base of the convection zone toward 
the surface. In Figure [6| we plot the maximum and the minimum of V(^^ in order to understand 
the differences in the asymmetries seen in Figure [Sj For each case, the maximum and minimum 
azimuthal velocities are measured over the whole domain but those extreme values are always 
reached inside or in the vicinity of the magnetic structures. Moreover, the various cases (except for 
Case Ob) only differ by the initial magnetic field strength of the loop and the differences in will 
thus be directly related to the differences in the magnetic properties. The most striking feature is 
the strong increase of the maximum velocity in all cases. This sharp increase happens in all cases 
when the loop reaches the mid-CZ, i.e., when the loop starts to get a significant curvature. If we 



look at Equation (19), we can understand that this sharp increase is due to the gradients with 
respect to the ^-coordinate of B^ and Bq which appear only when the loop becomes arched, i.e., 
after a certain rise time, contrary to the gradients of B^^ with respect to r and 0, which exist at 
the beginning of the simulation already. If we now turn to look at the minimum of v^^^ we realize 
that this change of slope is also visible at the same time but only for the strong field cases. They 
are indicated by the red arrows in panels (c) and (d). For the retrograde flow with respect to the 
rotating frame indeed, another ingredient plays a role, the conservation of angular momentum. As 
we said before, it will have a tendency to produce a negative which is dominant in the weak-field 
cases (panels (a) and (6)) but dominated by the Lorentz force terms in the strong field cases (panels 
(c) and (d)). As a consequence, strong gradients of the magnetic field and the associated magnetic 
tension will have the tendency to oppose the deformation of the loop while it rises and it explains 
why the asymmetry between the two legs of the loop is reduced when the initial field is strong 
enough (more than lO^G). This feature was already observed in thin- flux tube calculations (e.g.. 



Moreno-Insertis et al.|[T994{ Caligari et al. 1995) and a similar argument involving the stiffness 



of intense fields was used to explain the weaker asymmetry seen in strong field cases. We note 
moreover that since the rise time of strong tubes is smaller, the Coriolis force will have less time 
to act on the loop and the asymmetry will also be reduced because of shorter rise times. 

Going back to Figure |4} we see that when the loop forms and rises, some negative toroidal field 
is created (blue colors in the figures), while the initial toroidal field was purely positive. This can 
be explained by a local fi-effect: the azimuthal velocity gradient built at the periphery of the loop 
while it rises will shear the small poloidal field component present at the tube boundary mainly. 
This fi-effect, related to a shearing of poloidal field is one of the main ingredients of the dynamo 
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Fig. 6. — Evolution of the minimum and the maximum azimuthal velocity in the whole domain 
in the four cases, Case Ob (panel (a)), Case ITwP (panel (6)), Case 2TwP (panel (c)), and Case 
5TwP (panel (d)) of Figure [5| 



mechanism thought to be responsible for the generation of magnetic field in differentially rotating 



astrophysical objects (see for example Brunl 2004 ; Charbonneau 2005; Miesch et al. 2005 for a 



review on these processes in the solar context). The presence of such an effect could be surprising, 
especially in the untwisted case without convection, since no poloidal field nor azimuthal velocity 
gradients exist at the initial stage. However, while the loop forms, the magnetic field is no longer 
purely toroidal since the field lines get curved and will thus get a mainly radial component. At the 
same time, gradients of azimuthal velocity will be created as shown above. Both these effects are 
then favorable to get an fi-effect which will be able to create additional toroidal field of both signs. 



4.3. Structure of the Emerging Regions 

We now turn to investigate the characteristics of the individual bipolar regions emerging at 
the top of our computational domain in all situations. In particular, we focus on the tilt angle, the 
asymmetry in the field strength in the two spots and the morphological properties of the regions 
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which can be related to observations. 



4.3.1. Evolution of the Tilt Angle: Effects of Twist and Latitude 



We decide here to focus on four cases: Cases ITwP and ITwN and Cases 3TwP and 3TwN. 
The main difference between the two cases with right-handed twist is the field strength which 
is multiplied by 2 in Case 3 but buoyancy is kept the same through a change of initial entropy 
perturbation. Moreover, we know from previous studies (e.g., |Fan|[2008l ) that the sign of the initial 
twist will have a dramatic effect on the tilt angle of the emerging region. We address here this 
question by computing Cases 1 and 3 with a right-handed as well as left-handed twist. 
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Fig. 7. — Zoom on the emerging radial field at r = O.93i?0 for Cases 3TwP and 3TwN where 
So = lO^G. 

Figures [7| and |8] show a zoom on the emerging region in both cases by the time the apex of the 
loop reaches the top of our domain (O.96i?0). Blue indicates negative radial field and red positive. 
The first thing we note is the typical evolution of the inclination of the bipolar structure. Starting 
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from north/south because of the twist of the loop apex, it becomes more and more east /west as 
we see the legs of the loop emerging. In agreement with previous thin-flux tube calculations (e.g. 



D'Silva &; Choudhuri 1993), we find that a stronger initial field produces a smaller tilt angle since 
the magnetic tension will tend to oppose the tilting of the magnetic structure. This effect is visible 
when comparing the left panels of Figures [7| and [8] which differ by a factor 2 in initial field strength. 

The tilt angle of the bipolar structure and the effect of the initial twist is consistent with 
previous simulations, especially of Fan ( 2008| ). Indeed, in the left-handed case (which is thought 



to be dominant in the Northern Hemisphere), the tilt angle is not in agreement with observations 
which tend to show an inclination of the bipolar spots much closer to what we get in the right- 
handed twist. It was argued previously ( Fan||2008 ) that to reconcile simulations with observations. 



a very weak left-handed twist had to be introduced so that the Coriolis force acting on the loop 
would be strong enough to change the tilt angle from positive to negative sign. The problem is 
that such a weak twist, not sufficient to maintain a full coherence of the magnetic structure, leads 
to a strong flux loss and such structures would be likely to have difficulties to rise in a convective 
environment. 



If we look at Figure [s] where the initial field strength is decreased to 5 x lO^G but the amount 
of twist kept the same, we manage to get a case where the initial left-handed twist is strong enough 
to maintain the coherence of the loop and where the tilt angle of the emerging radial field gets 
in better agreement with observations. The Coriolis force has been able to act sufficiently on this 
loop to rotate the two opposite polarities clockwise and produce a final tilt angle of the same sign 
in both cases. 

To illustrate this. Figure [9] shows on the same plot the temporal evolution of the tilt angle 
measured in Cases ITwP and ITwN, at r = O.93i?0. The tilt angle was here determined by locating 
the peak values of each polarity and measuring the angle between the line linking them and the 
East- West direction. At the beginning of emergence, the angle is close to 90^ for both cases, but 
with opposite signs. As the emergence progresses, the absolute value of the tilt reduces in both 
cases. At this stage, the effect of the twist is dominant and corresponds to what is measured here. 
After about 14 days of evolution, the rapid change of tilt stops and the effect of the Coriolis force 
starts to be visible, especially on Case ITwN. As a consequence, the curve in Case ITwP seems 
to saturate around the value of 55^ and starts to increase again slightly because of the rotation of 
the spots, to reach a value of about 60^ after 26 days of evolution. In Case ITwN, the tilt angle 
continues to increase and crosses the ^ = line after about 15 days. At this instant, a change 
of sign has occurred and the tilt has thus become compatible with observations. A final value of 
about 15^ is reached after 26 days when diffusion then starts to act, which lies in the acceptable 
range of values when compared to observations at this latitude. 

We now turn to study the influence of latitude on the tilt angle. As we expect and as seen 



in Figure 10, the tilt angle depends on the initial latitude since the Coriolis force depends on the 



latitude. From the estimate of D'Silva &; Choudhuri (1993), the tilt angle is a decreasing function 
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Fig. 8. — Same as Figure [7]but for Cases ITwP and ITwN where Bq = 5 x lO'^G. 
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Fig. 9. — Measure of the tilt angle of emerging bipolar spots for Cases ITwP and ITwN, at 
r = O.93i?0. 
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of the cosine of the colatitude (or an increasing function of latitude), as seen in Figure 10 and 
as seen in observations. As we saw before in the left-handed case, the bipolar spot will have the 
tendency to rotate clockwise because of the Coriolis force acting differently on the two legs of the 
loop. As a consequence, the orientation of the bipolar spots in the case where it was introduced at 
high latitudes (60^) is close to north/south, and not only because we only see the twist of the field 
lines here but rather because the axis of the loop has been rotated clockwise by about 75^. When 
the initial latitude is 15^, the typical tilt angle we get is around 35^ while it increases to 60^ when 
a loop at initially 30^ is considered. Moreover, we note on this figure that the emergence occurs at 
slightly different locations in longitude. The center of the bipolar structure being located at around 
75^ in Case ITwPlS and 90^ in Case ITwPGO. This is due to the asymmetry of the loop which is 
much more pronounced in the low latitude cases since the action of the Coriolis force on both legs 
of the loop is stronger then. 

We note here that the values for the tilt angle may appear large compared to the classical 



observations (between 3 and 10^ typically) and to previous thin flux tube calculations (iD'Silva &; 



Choudhuri |1993[ [Caligari et al. |1995| ) but as shown before (on Figure [9]), an initial left-handed 



twist, associated with a slow rise allowing the Coriolis force to act on the loop, is able to produce 
tilt angles of about 15^ in a case where the initial latitude was 30^. This, in turn, is very compatible 
with observations at the photosphere. It is possible here to make such a comparison if we make the 
reasonable assumption that the Coriolis force will not have the time to modify the loop orientation 
and thus the tilt angle during the fast rise through the last few percents of the convection zone. 

We shall keep in mind those results for the convective cases where the spherical shell will now 
be differentially rotating. Since the tilt angle is also a function of the local rotation rate, the tilt 
angles we will get in the convective cases may well be different. 



4-3.2. Morphology of Emerging Regions: Tongues and Necklaces 



If we concentrate again on Figure [7| we immediately notice the typical "tongue" -shape of 
the emerging region, especially at a deeper radius r = O.87i?0. They are indicated by red and 
blue arrows in the bottom left panel of Figure [?[ Those elongated structures, also seen in the 



observations (Lopez Fuentes et al. 2000), in Cartesian simulations (Fan 2001; Archontis &; Hood 



2010|) and studied in details in |Luoni et al.| ( |2011| ), are mainly due to the twist of the field lines 
showing in the map of the radial magnetic field. In particular, the orientation of the tongues 
(or the direction in which they extend) give an indication of the sign of the twist. In the lower 
left panel of Figure [7} we thus see that both positive and negative polarities extend on their left 
side when it is the opposite for the lower right panel where the sign of the twist has changed. 
Another interesting feature here is the appearance of sharp structures around the emerging region, 
forming an annular shape that we will from now on call a "magnetic necklace". This annular 
shape of emerging active regions is also often seen in observations but has never been reported as 
identified magnetic structures, different from the so-called tongues. This is true for observations 
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Fig. 10. — Measure of the tilt angles of cases where the loop was introduced at the latitudes 
of 15^ (Case lTwP15), 30^ (Case ITwP) and 60^ (Case ITwPGO). The Northern Hemisphere is 
represented centered at 90^ in longitude and 45^ in latitude. 



with experiments and instruments as varied as the Flare Genesis Experiment ( Pariat et al.||2004 ), 
Hinode dOtsuji et al.|[20TT1 ) or SOHO/MDl ( |Liu fc Zhani|[2006l ) . 
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Those structures are clearly visible in the lower panels of Figure [7| but do not really show in 
Figure|8]because of the low field strength and the rather large magnetic diffusivity in our simulation. 
Indeed, since those structures are organized at rather small scales, it is necessary to have low enough 
diffusivities to be able to see those structures at the time of emergence and during their rise. As 
a consequence, we decided to compute Case ITwN with a lower magnetic diffusivity (we choose 
Pm = 5 here) to see if these small-scale structures could be recovered in this case. Indeed, as shown 
on Figure pT| the magnetic necklace is visible on the zoom in the (6, (f)) plane on the emerging region 
at O.93i?0 and is also seen in a meridian cut of the loop during its rise (see red and blue arrows in 
the left panels). They consist in sharp elongated structures (not to be confused with the tongues) 
lying at the periphery of each polarity and of the same sign that create a ring inside which the 
two main regions of strong radial field concentrate. If we compare the regions of appearance of 
these sharp structures around the bipolar spot with the norm of the vorticity (or enstrophy), it is 
rather clear that those structures correspond to a velocity shear which produces a strong vorticity. 
As shown and studied in Emonet & Moreno- Insertis (1998), a sharp interface is indeed formed at 
the boundaries of the flux rope while it rises. This region is the site of vorticity generation via the 
presence of a magnetic field during the whole loop evolution. Another source of vorticity exists 
within the tube interior and is linked to the gravitational torque applied between the tube center 
and its periphery. This process creates the counter-rotating vortices which can completely break up 
the tube during its rise if the twist (and thus the magnetic tension) is not strong enough, as seen in 
Section [4. Looking at the four panels of Figure [TTJ we understand that the necklace, associated 
with vorticity generation at the loop edges, will be visible only after a significant part of the loop 
has emerged. Indeed, since the main vorticity generation is located at the sides of the loop, we need 
to wait until the axis of the loop crosses the surface r — O.93i?0 before being able to observe the 
sharp magnetic structures. By that time, the tongues have thus already significantly faded away 
and reconcentrated around the two polarities. We can thus expect to see those structures in actual 
observations only after a certain amount of time after the first signs of emergence. 



4-3.3. Field Strengths and Fluxes 

We end this section by noting another difference between Cases 1 and 3. In Case 3, the loop 
seems to maintain more coherence while it rises, both polarities are well-separated, even at the end 
of the simulation when the bipolar region reaches its maximum extension. This is less the case for 
the weaker field of Case 1. Indeed, both polarities seem to get mixed and the radial magnetic field 
lines get more easily advected by the velocity field created by the loop itself (we remind that there 
is no convection in this case). At the end of the simulation, the magnetic field gets advected at the 
boundary of the positive radial velocity pattern it creates and finally diffuses away. 

If we now focus on this more realistic Case 1 (with a field strength of 5 x IC^G, and a typical 
rise time of more than 10 days), we can move to a more quantitative study of the evolution of 
the intensity of the radial field emerging in each spot and of the amount of flux it represents. 
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Fig. 11. — Radial magnetic field and norm of the vorticity in a (0, plane (upper panels) and a 
(r, 0) plane (lower panels), focusing on the sharp structures (indicated by arrows) which surround 
the emerging regions. 



The results of this study are shown in Figure 12, Panels (a) and (6) show the evolution, during 
emergence at r = O.93i?0 of the radial field and radial flux in Case ITwP and panels (c) and (d) in 
Case ITwN. In both cases at the beginning of emergence, the negative field in the leading polarity 
dominates over the trailing polarity (by about 35%), while the amount of flux is the same in both 
polarities (since all emerging field lines of one polarity connect to the other one in the domain of 
integration). This suggests that the leading polarity is stronger in intensity and more concentrated 
and the trailing polarity is weaker and less coherent, which is in agreement with observations of 
emerging active regions (e.g., [Kosovichev fc Stenfio 2008). Moreover, this feature was also first 
seen in Cartesian simulations and explained in Fan et al. (1993j) by the evacuation of plasma out of 
the leading side of the loop into the following side. The leading side then shrinks and the magnetic 
pressure increases to balance the same external pressure, resulting in a more intense field in the 



less extended leading polarity. As Fan et al. (1993) points and as discussed also in Caligari et al. 



(1998), the anchoring of the field lines outside the buoyant part of the loop is a key ingredient 
to produce a strong asymmetry in field strength. It is also what we find here in Case ITwN, for 
example, where a strong asymmetry in inclination develops because of the anchoring of the field 
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Fig. 12. — Radial field and unsigned radial flux for Case 1 with right-handed (panels (a) and (6)) 
and left-handed (panels (c) and (d)) initial twist. The values of the field are measured at r = O.93i?0 
in the center of each polarity where the field is strongest and the fluxes are the integrals of the 
radial field contained in the trailing (or leading) polarity through a surface enclosing the emerging 
region. 



lines outside the most buoyant part and where we find a leading leg about 1.4 times more intense 



than the trailing leg on average (see panel (c) of Figure 12 



We note that the values of the flux are quite strong compared to observations (which report 
values of about 10^^-10^^ Mx). This is due to the scales of our magnetic structures which over- 
estimate the size of a flux tube initially located at the base of the convection zone. The initial 
flux in the magnetic structure in Case 1 is indeed of 6 x lO^^Mx. Simulations with flux tubes of 
radii about 10 times smaller, numerically very challenging, would lead to values for the fluxes much 
closer to observations ( Schrijver fc Harvey||1994 ). It is however important to note that despite the 
discrepancies between the absolute values of the fluxes, due to numerical limitations, the shape 
of the temporal evolution of the fluxes is in good agreement with observations. Indeed, the flux 



in each polarity of the bipolar active regions is observed on MDI magnetograms for example (van 



Driel-Gesztelyi et al.||2003D to increase sharply in the first days of emergence, before saturating and 
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slowly starting its diffusive decay. We get exactly the same kind of evolution in our simulations, 



visible especially in panel (b) of Figure 12 The flux strongly increases from to 1.5 x lO^^Mx in 
the first 5 days before saturating around this value. In the figure and especially in panel (d), we 
distinguish first a slow rising phase between 8 and 9 days followed by the very sharp increase of the 
amplitude of the flux after 10 days of evolution. The first phase is related to the very beginning 
of the emergence when only the apex of the loop is visible and the only contributor to the flux 
is the radial magnetic field coming from the twist of the field lines. As the emergence proceeds, 
however, the feet of the loops where the strong radial magnetic field lies start to emerge and form 
the bipolar region. It is thus this fast emergence phase of the arched structure which produces 
the sharp increase of the simulated (and possibly of the observed) magnetic fluxes. As far as the 
magnetic field strength is concerned, we can wonder what typical horizontal field strength we would 
get at the photosphere, starting with our flux tubes located at the base of the CZ, in the same way 
as Cheung et al. (2010). If we follow the evolution of the axial field intensity in time (and thus in 
radius), we find that the scaling with the reference density is close to B (x {p-\r p) (similar to what is 
found in |Pinto fc Brun|2012 ), corresponding to an expansion of the loop in the directions transverse 



to its axis. At the top of our domain, the typical toroidal (or horizontal) field is of the order lOkG 
(for Case ITwP for example) and the density p of the order 4 x 10~^g.cm~^. As a consequence, 
if we consider a density of 4 x 10~^g.cm~^ at the photosphere, the horizontal field strength would 
reduce to 1 G, much weaker than the typical observed values (around lOOG). However, the expan- 
sion in the last few percent of the convection zone might be modified (becoming mostly horizontal, 
for example, due to a decrease of the plasma j3) and then the scaling with density could change 



drastically and become B oc ^/p as shown in Cheung et al. (2010). We should then keep in mind 



that we are not modeling the uppermost layers of the convection zone and that it is not relevant 
here to relate directly the field strengths in our bipolar regions to the field strengths which would 
be reached at the photosphere. 

Comparing the left-handed and right-handed twist cases, we find that the evolutions, especially 
for the radial field, are slightly different. Indeed, the dominant leading polarity at the beginning 
of emergence in Case ITwP (panel (a)) saturates at the value of 2700G and then drops to the 
value of 2300G while the trailing polarity continues to increase in intensity. At t = 14 days, the 
positive trailing polarity thus becomes dominant, before decreasing at t = 21 days. This is not the 
case for Case ITwN (panel (c)), where the leading polarity always stay dominant in strength and 
approximately equal to the trailing polarity in flux. This could be due to the asymmetry built up 
during rise, which is smaller in the left-handed twist case. This is indeed seen on cuts at constant 
latitude as in Figure [s] (not shown here). As a consequence, in Case ITwP, the asymmetry between 
the legs is strong, the leading polarity reaches its maximum strength and thus starts to diffuse 



away before the trailing polarity. This explains the decrease in amplitude in Figure 12 first seen 
for the negative field and then for the positive field. In Case ITwN, the peak positive and negative 
radial fields are measured at r = O.93i?0 at approximately the same time since the asymmetry is 
less pronounced and they both saturate at the same time (at the value of 1900G for the positive 
polarity and 2500G for the negative one), before starting their diffusion. 
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5. EVOLUTION OF A SINGLE LOOP WITH CONVECTION 

We now turn to investigate the behavior of a loop similar to the cases presented above, but 



inside a convective background. As in the global flux tubes simulations (Jouve &; Brun 2009), the 



buoyancy force will compete with the downflows to drive the loop towards the top of the computa- 
tional domain and the emerging region will have significantly different characteristics than in the 



isentropic cases. Figure [13] shows the typical three-dimensional evolution of an fi-loop evolving in 
a convective environment, this corresponds to Case 2TwP. We only see here the evolution of the 
toroidal field but a turbulent convective velocity field exists all around it. Moreover, differential 



rotation and meridional circulation have well-established profiles here, as stated in Section 2.3 



Time=4.1 days 



Time=6.2 days 



Time=83 days ^ Time=10.8 days 



Fig. 13. — Volume rendering of while the loop rises through a convective layer. As in Figure 
[4} the loop is viewed from below, looking up at the North Pole. The deformation of the loop by 
convective motions is clearly visible here when compared to Figure |4| 

We recover some aspects of the cases studied in Section |4j like the asymmetry between the 
legs and the creation of a negative toroidal field at various regions at the loop periphery. Again, a 
local fi-effect is at work here, with a slightly more complex structure since the velocity field is now 
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organized at different spatial and temporal scales. But the main ingredients (azimuthal velocity 
gradients and radial and latitudinal components of the magnetic field) are still the same and able 
to locally create negative toroidal field. 



5.1. In the Bulk of the Convection Zone 
Case ITwP CaselTwP 




Fig. 14. — Cut at the latitude of 30^ of the magnetic energy (contours) superimposed to the 
background radial velocity (blue meaning downflows and red meaning upflows) for Cases ITwP 
and 2TwP. The color bar indicates the values taken by the radial velocity. 



In such fully convective setting, we expect the initial field strength in the flux tube to be an 
important parameter for the evolution and rise of the magnetic field. Indeed, as shown in previous 
work involving a uniformly buoyant flux tube ( Fan et a"L||2003 ; Jouve fc Brun||2009 ) , the downflows 
and upflows of the convective layer control the rise velocity of the tube and are sometimes even able 
to pin it within the convection zone, resulting in structures emerging at specific longitudes. Here 
again, we get the same kind of competition between convection and magnetic buoyancy, which can 



be seen on Figure [Mj This figure shows the evolution of the loops of Case ITwP and Case 2TwP. 
They correspond to the second and third panels of Figure [5j for which the asymmetry was clearly 
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modified when the parameters were changed such as to have more intense field. In Case 1, the 
initial field strength is 5 x lO^G, i.e., slightly below the equipartition field strength of 6 x lO^G, 
corresponding to the magnetic field in equipartition with the strongest downdrafts at the base 
of the CZ. In Case 2, the initial magnetic field strength is lO^G, corresponding to 1.7 times the 
equipartition field. We superimpose here the magnetic energy contours to the background radial 
velocity (representing the convective motions). For Figure [5j we chose to represent the loops at 
about t = 10 days since their position in the bulk of the CZ was not very different for the isentropic 
cases. Here, at about t = 10 days, the loops have not reached the same height, the apex of the 
loop of Case 2TwP has reached r = O.93i?0 after about 8 days already, whereas the loop of Case 
ITwP reaches the same height at t = 13 days and was thus significantly slowed down compared 
to the isentropic case. This can be explained by the structure of the flow in which the loops are 
embedded. In Case 2TwP and as seen especially on the panel at time t = 8.11 days, the velocity 
field associated to the Lorentz force is dominant. At the location of maximum magnetic field, a 
strong upflow thus develops (a maximum positive radial velocity of 350m. s~^ is reached here), the 
loop is dragged upward and the influence of the downflow shown by the black arrow is limited. 
The rise of the loop is thus widely promoted by the strong upward convective cell produced by the 
intense magnetic field. On the contrary, in Case ITwP where the magnetic energy content (and 
thus the buoyancy) is weaker, the radial velocity created by the presence of a magnetic field is of the 
same order as the background convective motions and the downflows (as the one indicated by the 
arrow) are much more efficient at maintaining the magnetic loop under the surface and preventing 
the axis from emerging. 

Another difference we note here compared to the isentropic cases is the degree of asymmetry 
between the trailing and leading legs. The strong asymmetry visible in Case ITwP when the loop 
was embedded in an isentropic background is still present at t = 13.3 days but less obvious since 
the advection by the convective motions has also mixed the field lines and significantly modified 
the structure. However, in Case 2TwP at t = 10.4 days, the same degree of asymmetry as in the 
isentropic case is found, showing again that this regime is dominated by the Lorentz force, contrary 
to Case ITwP where the magnetic structure is essentially frozen in (advected by) the convective 
motions. 

In order to gain a better insight on why we get such a different behavior between what we 
will now call the convection dominated regime of Case ITwP and the magnetic field dominated 
regime of the other cases, we plot the evolution of the density deficit inside the loop at its apex as 



a function of time in Cases 1 and 3. The results are shown in Figure 15, They correspond to Cases 
ITwN and 3TwN but the results are similar for the right-handed twist cases. Our results here 
and especially the distinction between convection and magnetic field dominated cases are similar to 



what is obtained in the recent thin-flux tube simulations of Weber et al. (2011). Indeed, they find 
that an initial field strength above lO^G ensures an evolution dominated by magnetic buoyancy 
while a loop initialized with Bq = 1.5 x lO^G is strongly influenced by convective motions. Similar 
values are found here. 
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Fig. 15. — Density deficit at the apex of the loop during its rise in Cases ITwN (dotted hne) and 
3TwN (sohd hne). The internal density is measured inside the tube apex, where the magnetic field 
is maximum and compared to the external density which is taken to be the average value of the 
density at the corresponding depth. 

F igur e [15] shows the evolution of the density deficit inside the loop compared to its surroundings, 
for Case 1 (dotted line) and Case 3 (plain line). This density deficit represents the buoyancy of the 
structure and thus its capacity to reach the top of our computational domain. Both loops started 
with approximately the same density deficit, as seen at the beginning of the evolution. The two 
cases first follow the same amount of buoyancy loss during about half a day and then the two 
cases start to have different behaviors. Indeed, in Case 1, the loops continues to smoothly lose its 
density deficit while the simulation proceeds and follows the exponential decay of a typical solution 
of diffusion equation. After t = 10 days, the density fluctuations in the bulk of the convection zone 
are of the same order as the density deficit in the loop and the magnetic structure stops being 
buoyant. We also note that convective motions act on the magnetic structure to create smaller 
scale regions which will thus diffuse faster and lead to a smaller field strength and thus smaller 
buoyancy. At t = 15 days, the loop has still not reached the top of the domain and the only reason 
why we will get some emerging magnetic field in this case is that upflows will drag the loop upward. 
The coherent bipolar structure will then be mostly lost when the magnetic structures will emerge 
at r = O.93i?0, as we will see in the following section. In Case 3 on the contrary, the evolution of 
the density deficit is quite different. The major trend is of course a decrease of the buoyancy of 
the structure but this decrease does not follow a decaying exponential as we had in the previous 
case. Here, strong currents generated by the higher field strength (initially lO^G) will play a role 
in the energy balance of the loop and the diffusion of entropy will be less efficient in this case. As 
we see on the figure, the buoyancy of the loop of Case 3 always stay higher than in Case 1. As a 
consequence, at t = 15 days, a significant part of the loop has emerged at the surface as a coherent 
bipolar structure. The loop thus keeps both its coherence and its buoyancy for a longer time and 
the structure of the emerging regions will tend to be less influenced by the convective motions. 
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5.2. Structure of the Emerging Regions 

As stated in the introduction, one of the main questions we wish to answer in this work is 
the fohowing: "what is the influence of convective motions and large-scale flows on the structure 
of our emerging magnetic field?". Although our top boundary condition is still at 28Mm below 
the surface and that density and pressure drops may well cause large modifications of the observed 
magnetic field, we can get some insight on what type of active regions may be created at the solar 
surface from our numerical simulations. 



5.2.1. Convection versus Magnetic Field Dominated Regimes 

In this section, we focus on the convection dominated regime of Case ITwP (and ITwN), 
where the emerging magnetic field will be strongly affected by the convective motions. Figure [16] 
shows four snapshots of the emerging region zoomed in latitude and longitude, at r = O.93i?0. The 
radial velocity (colors) represents the convective motions (blue downdrafts and yellow upflows) 
and we superimpose the emerging radial field (red contours indicating positive polarity). At the 
beginning of emergence at about t — Q days, a well-defined bipolar structure can be identified, 
at the latitude of 30^, indicating that the rise was indeed radial and not parallel to the rotation 
axis. We note already on this snapshot that the strong downflow located at the latitude of 30^ 
and the longitude of 105^ strongly interacts with the negative polarity and is responsible for its 
cashew-nut-like shape. As the emergence proceeds, the well-identified bipolar structure is replaced 
by a much more complex field structure. The negative polarity gets squeezed into the previous 
strong downflow, while regions of positive polarity start to emerge all around. The concentration 
of magnetic field in the downflow lanes is even more visible on the third snapshot dX t — 16.1 days 
where the field concentration seen on the previous panel continues to accumulate in the downdraft 
while a new emerging bipolar structure appears at the longitude of 90^ where a large upflow was 
already present and responsible for this new emergence. We then note that this new region gets 
quickly advected toward the boundary of this large upflow to, in turn, get concentrated into the 
neighboring downdrafts, leading to an evolution very close to what is observed in simulations of 



magnetoconvection ( , Weiss et al. 1996). In this case, convective motions were thus very weakly 



affected by the emerging magnetic field and completely dominates its evolution. 

The evolution shown in the previous figure is quite different when cases with stronger buoyancy 
and a stronger initial field strength are considered. We now turn to investigate the structure of 
the magnetic field emerging in Cases 2TwP and 4TwP which initially possess the same buoyancy 
but different field strengths (lO^G for Case 2TwP and 1.35 x lO^G for Case 4TwP). These values 
respectively correspond to 1.7 and 2.3 times the equipartition field strength. The results are shown 
in Figure [T7| Four snapshots are shown, two for each case, at approximately the same times, t — Q 
days and t — lldays. Case 2TwP is interesting since it is an intermediate case between convection 
dominated and magnetic field dominated evolution. Indeed, at the beginning of emergence, a bipolar 
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Case ITwP: Convection-dominated 




Fig. 16. — Zoom at r = O.93i?0 on the emerging radial field (red contours represent positive B^) 
inside the convective environment (yellow colours represent upflows), in Case ITwP. The typical 
magnetic field strength in these snapshots is about 2000G. 



structure emerges at about 95^ as expected (and as found in the isentropic cases) but another region, 
which is in fact another part of the loop, also emerges at the longitude of 105^. This is due to the 
strong downflow lane at about 100^ of longitude and its neighboring upflow which will have the 
effect of dragging the loop upwards in the upflow and pinning it down in the downflow (those flows 
are also visible in Figure 14). As a consequence, two separate regions seem to emerge when they 
are part of the same fi-loop but strongly influenced by convection. As the simulation evolves and 
emergence proceeds, we recover that the magnetic field gets advected in the downflow lanes at the 
largest longitudes but shapes as a typical active region at around 90^. Indeed, the tongue-shape 
appears here, although modified by the convective motions. This typical case, intermediate between 
a convection dominated and magnetic field dominated evolution, is also interesting since it directly 
relates to observations of emerging regions. Indeed, the regions of strong magnetic field, in the 
center of the active regions, are only marginally affected by the convective motions, as opposed to 
the surrounding plage regions whose structure is affected by the supergranular scales of convection. 
We could thus consider that the parameters of Case 2 enable us to simulate an emerging region with 
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characteristics typical of an observed active region, in particular its field strength in each polarity, 
its competition with convective motions, its tongues or its tilt angle. 

If we now look at Case 4TwP where the initial magnetic field was stronger, the emergence is 
much more similar to what we had in the isentropic cases. The orientation of the unique (unlike 
Case 2) structure is first North-South to become more and more East- West, the tongues are easily 
identified, a strong upflow is created by the loop itself and dominates its evolution during the 
emergence process. We observe in Figure [T7| that the apex of the loop emerges at the longitude 
of 95^ and that no radial field is visible at 105^ in Case 4, contrary to Case 2. This is due to the 
fact that in Case 2, the advection of the magnetic structure by the convective motions significantly 
modifies the shape of the loop and thus the asymmetry between the trailing and leading legs. In 
Case 4, the evolution is again closer to the isentropic case and convection has little influence on 
the loop and is thus not strong enough to produce an emerging region within the upflow located at 
105^. To conclude on these cases, although Case 2 is more influenced by convection, a clear impact 
of the emerging structure on the velocity field is visible in both cases, especially at t = 11 days 
where a strong upflow is created at the very location of emergence. 
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5.2.2. Effect of Twist and Latitude 



In this last section, we investigate the effect of changing the sign of the twist of the initial 
magnetic loop, as well as modifying its latitude of introduction. We saw in the isentropic cases 
that the sign of the twist would influence the tilt angle of the emerging region but that a tilt 
compatible with observations could still be found with an initial left-handed twist. We now address 
the question of the possibility to get tilt angles even more compatible with observations when the 
magnetic loop is embedded in a convective environment. 
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Fig. 18. — Temporal evolution of the radial magnetic field intensity in the trailing and leading 
polarities at and right after emergence, for cases ITwP (panel (a)), ITwN (panel (6)), 3TwP 
(panel (c)) and 3TwN (panel (d)). 



Before specifically studying the tilt induced by the twisted structures, we investigate the impact 
of a negative initial twist on the intensity of the emerging radial field. Indeed, it was already found 
in the isentropic cases that the strength of the emerging polarities and the positive and negative 
fluxes had a different evolution when a left-handed or right-handed twist was considered. It is still 



the case when the loops are introduced in a convective shell. Figure [T8| presents the evolution of the 
radial field of the leading and trailing polarities at the time of and after emergence for Cases ITwP, 
ITwN, 3TwP and 3TwN. Similar pictures for Cases ITwP and ITwN were shown in Figure [12 
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for the isentropic simulations. We can thus first compare the convective and isentropic situations 
for Case 1. The first thing we note and which is true both for the right-handed and left-handed 
twist, is that the loop in the isentropic case had a quicker rise than in the convective case. Similar 



results were found for uniformly buoyant flux tubes in Jouve & Brun (2009). Indeed, for example 



at t = 10 days for the isentropic cases, the typical field strength already reached 2000G whereas it 



barely gets to 500G in the convective cases. This is related to the explanations of Section [5T] where 
the buoyancy loss was shown and thought to be partially due to the convective motions creating 
smaller scales and thus stronger dissipation. 

Another striking feature is the strong increase of negative Br in Case ITwP (panel (a)) at 
t = 13 days and of positive in Case ITwN (panel (6)) at about the same time. This is mainly 
due to the fact that these cases are dominated by convective motions. Indeed, as we saw in FigurjTHl 
the negative polarity in Case ITwP between t — 12 and t = 16 days gets strongly squeezed by 
the convective motions in the narrow downflow lanes. This compression of magnetic field tends 
to produce stronger concentration of fields, explaining the strong increase of negative B^ in Case 
ITwP and then its sharp decrease due to fast dissipation of small scales. A similar situation 
happens in Case ITwN for the opposite dominating polarity. In Cases 3TwP and 3TwN (panels 
(c) and (g?)) where the initial magnetic field is twice as large, such sharp increase and decrease 
of the two polarities are not visible since the loop emerges as a global magnetic structure and 
is much less influenced by convective motions. It is worth noticing that we may have a possible 
way here to explain the formation of complex active regions from where powerful (X-class) flares 
could originate. Indeed, active regions with a complex topology could be created by the turbulent 
convective motions acting on moderately strong flux tubes (of approximately 5 x lO^G at the base 
of the CZ). By the time they reach the top of the convection zone, the flux tubes have been strongly 
deformed and the subsequent emerging radial field consists in several patches of opposite polarities, 
as seen in Figure [iGj We could argue that this complex active region, whose progenitor is a relatively 
weak flux tube, would not be able to produce very energetic flares. However, when panel (a) of 



Figure [T8| is considered, we conclude that the converging motions at the sub-photospheric level are 
able to re-concentrate the magnetic structure in the strong downdrafts, in order to locally produce 
a strong increase of the radial field. In particular, we note that the maximum positive radial field 
is about twice as large as the negative S^, which was not the case in the isentropic case where 
the re-concentration by convective motions could not occur. This scenario may thus explain the 
injection of energy in the corona by convective motions acting on complex active regions originating 
from initially weak flux tubes. 

As far as the tilt angle is concerned in these simulations, the influence of the initial twist, 
which is not included in the thin-flux tube calculations, is again crucial. If we follow the position 
in latitude of the positive trailing and negative leading polarities (not shown here), we find that, at 
the beginning of emergence, the leading polarity is located at higher latitude than the leading one, 
unlike in solar observations. However, in Case ITwN, the clockwise rotation due to the Coriolis 
force acting on the legs of the loop tends to reverse this situation about 5 days after emergence. 
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implying a leading polarity at slightly lower latitude than the trailing one and thus to a tilt of 
about 5^, in agreement with Joy's law. This is similar to what was found in the isentropic Case 
ITwN. Case 3TwN does not follow the same evolution. Indeed, the orientation at emergence is 
kept even after 10 days after emergence before the magnetic field is swept towards the downflow 
lanes and dissipates. Again, we need a loop that is sufficiently influenced by the Coriolis force to 
get a tilt in agreement with observations with an initial left-handed twist believed to dominate in 
the Northern Hemisphere. 

Of course, not only will the initial twist play a role in the tilt of the emerging region but also 
the latitude of introduction and emergence, as already pointed out in the isentropic cases. We thus 
computed extra cases as in the isentropic simulations. Cases 2TwP15 and 2TwP60. We did not 
choose the same initial parameters as in the isentropic cases to study the effect of latitude, since 
Case 1 was computed in the isentropic case and it is believed here to be influenced too much by 
convection to produce a significant trend and thus to allow drawing any conclusions on the tilt at 
the time of emergence at O.93i?0. Case 2TwP is the most interesting since it is still influenced 
by convection (as we saw in the previous section) without completely loosing its characteristic 
bipolar structure when the loop emerges at the top of our domain. The results of these extra 
calculations are shown in Figure [19] Both the emerging radial magnetic field and a measure of 



the tilt angle are represented, similar to Figure 10 We argue that some effects of the differential 



rotation are visible here. First, the rise velocities are different from one case to another, like the 
isentropic simulations but the differences are even more pronounced here. To be more precise, the 
characteristic rise time for the loop initially located at a latitude of 60^ is around 8 days, when it 
is of more than 11 days for the case at latitude = 15^. The figures represent the emerging region 
about 1 to 2 days after the first signs of emergence at r = O.93i?0. This slow down of the loops 
located at lower latitudes is related to the rotation rate which is stronger at lower latitudes. Indeed, 
as shown already in |Jouve fc Brun ( |2007 ), the modified buoyancy force in the presence of rotation is 



proportional to Ap(^ — r sin^ 9fl'^) and the rise velocity of the loop will thus be less at low latitudes 
and high rotation rates, which is typically the situation for Case 2TwP15 represented here in the 



first panel of Figure [T9| A second effect of the differential rotation which can be noted concerns 
the longitude of emergence. In these differentially rotating cases, the loop which was introduced 
at higher latitudes has drifted in longitude by about 30^ (it was introduced at 100^ and emerges 
at 70^) when the loop in Case2TwP15 has only drifted by 5^. This is quite different from what 
we had in the isentropic cases where the drift in longitude was more important for the case at low 
latitudes. This was due to the stronger asymmetry between the two legs of the loop which appeared 
in the low latitude case. Here, the different drift values could well be due again to the differential 
rotation. Indeed, since the higher (lower) latitudes rotate retrograde (prograde) compared to the 
rotating frame, it is not surprising that the loop at high latitudes will be globally shifted at lower 
longitudes compared to the low latitude cases. This is thus a clear effect here of the large-scale 
differential rotation acting on our rising fi-loops. Finally, the idea here was to compare the tilt 
angles we could get in these convective cases and compare those values to observations. It is more 
complicated here to identify well-defined concentrations of opposite polarities, since the convective 



-39- 



motions have influenced the loops during their rise to sometimes produce several concentrations of 
radial magnetic field of both polarities. However, if we consider that the tilt angle is defined as the 
angle between the East- West direction and the line joining the maximum of both polarities, we are 
still able to measure a tilt in those 3 cases. We find a value of about 20^ for Case2TwP15, 40^ for 
Case2TwP and 45^ for Case2TwP60. The trend is thus still to increase with increasing latitude 
but it is less obvious than in the isentropic cases. We note here that the tilt angle distribution in 
the simulations of Weber et al. (2011) was also much more spread in the convective cases than in 
the isentropic ones, which is consistent with our simulations. This is an indication here that it is 
not only the rotation and the associated Coriolis force acting on the loop which will determine the 
tilt angle observed at the solar surface but horizontal buffeting by convection could also strongly 
modify this measurement, even if the trend compatible with Joy's law is still found in those cases. 
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Latitude=15^ 




Time=12.5 d 

Latitude=30" 




Time=11.07 d 

Latitude=60'' 




Time=10.2 d 



Fig. 19. — Same as Figure 



but for the convective Cases 2TwP15, 2TwP, and 2TwP60. 
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6. Discussion and conclusions 

In this work, we presented the dynamical evolution of twisted fi-loops embedded into a convec- 
tively stable and unstable unmagnetized environment, to mimic the behavior of a strong toroidal 
magnetic field generated by dynamo action at the base of the solar convection zone and rising 
toward the photosphere. Those simulations allow to take into account the full effects of sphericity, 
mean flows and convective motions on the dynamical evolution of a flux rope with a finite radius 
and thus contains much more physics than previous thin flux tube calculations. The first part of 
this work was dedicated to establish the important parameters to get a coherent radial rise in the 
isentropic case. We find that a loop with an initial magnetic field of 5 x 10^ G and with an initial 
twist of less than 1 turn along the loop extension (15^ here) is able to rise and emerge at the top of 
our domain to create bipolar regions. Asymmetries develop between the leading and trailing legs of 
the loop due to the combined effect of angular momentum conservation and Lorentz forces and it 
is found that this asymmetry is strongly reduced when magnetic structures as strong as 1.5 x lO^G 
are considered. The emerging regions obtained in these isentropic cases have several properties 
in agreement with observations. In particular, well-identified tongues develop during emergence. 



mainly due to the initial twist of the magnetic field lines, also seen in observations (e.g. Lopez 



Fuentes et al. 2000). For the first time, we point the appearance during flux emergence of sharp 
ring-shaped magnetic structures around the bipolar region, that we call "magnetic necklaces" and 
which are due to the creation of vorticity at the loop periphery. They are also very often seen in 



observations of large-scale flux emergence (especially in Liu fc Zhang||2006 with MDI data) and are 



very similar in shape to the ones we get in our simulated active regions. Moreover, we show here 
that it is possible to get a case where the initial left-handed twist is strong enough to maintain the 
coherence of the loop and where the tilt angle of the emerging radial field agrees with Joy's law (in 
contrast to what is found in Fan||2008 ). To do so, the rise time needs to be long enough so that the 



Coriolis force acts sufficiently on the loop to rotate the two opposite polarities clockwise and then 
produce a final tilt angle of the correct sign. 

We note that if much thinner flux tubes were considered, the various parameters quoted here 
to get bipolar regions with properties in agreement with observations would probably have to be 
revised. The results do not change significantly if a tube four times smaller is introduced as long 
as the magnetic diffusivity is also reduced to keep a rise time small compared to the diffusion time 
of the magnetic structure (see Jouve fc Brun||2009| for a discussion on the effect of diffusion on this 



type of simulations). However, if much thinner tubes are introduced (say 10 to 20 times smaller 
as can be expected in the real Sun), it is likely that the results may differ. In particular, the drag 
force coming from the pressure difference between the up and down stream sides of the turbulent 
wake behind the rising loop may be much bigger when the loop radius is reduced. This was shown 
in Batchelor (1967) for incompressible flows passing a rigid cylinder under large Reynolds number 



conditions and shown to be still true for two-dimensional simulations of buoyant loops (Emonet 



&; Moreno-Insertis 1998). This will in turn modify the rise velocity and the competition between 



magnetic fields and convective motions. For example, a smaller rise velocity leaves more time for 
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the Coriolis force to act but at the same time, an increased drag directly opposes the tilting motion, 



resulting in a net decrease of the tilt angle of loops of smaller radii, as shown by D'Silva fc Choudhuri 



(1993). This work would then benefit from studying the effect of loops of smaller radii but this 



question was not addressed here. Pinto & Brun (2012) also find that the wake structure is modified 



by the presence of a background dynamo field, possibly due to continuous field reconnection as the 
tube rises. 

In the convective cases, the rise velocity and characteristics of emerging regions are strongly 
affected by the convective motions when loops of less than lO^G are considered. However, emerging 
regions with the correct orientation and a dominant leading polarity are still found in these simu- 
lations where the full effects of convection, rotation and sphericity are taken into account. Making 
direct comparison with observed photospheric field intensity in active regions is not straightfor- 
ward since the scaling of the magnetic field strength with density (which is linear in our cases) may 
change drastically when reaching the uppermost layers of the convection zone dCheung et aL||20Tol ). 
However, local helioseismology techniques seem to start to provide us with constraints on the field 
intensity deeper down (|Ilonidis et al.||2011 ; Braun||2012 ; Ilonidis et al.||2012 ), much closer to the top 
of our computational domain. We may well soon be able to relate much more directly the results 
of our three-dimensional simulations to actual observations. 

If the origin of the tilt angle of active regions has not been completely related to one particular 
property of the emerging loops yet, it is however undeniable that the initial twist of the field 
lines is an important ingredient. It is thus necessary to understand the generation not only of 
buoyant magnetic field at the base of the convection zone but of twisted buoyant structures. It 
is believed that such twisted loops could be created by the advection of poloidal field lines by 



the initially rising strong toroidal field subject to magnetic buoyancy instabilities (Longcope et 



al. 


1996 


; Favier et al. 


2012) 



buoyant, gaining some twist and being able to rise coherently to the top of the convection zone are 



in progress (Nelson et al. 2011, 2012). We moreover note that such simulations are able to mimic 
the dynamical evolution of a magnetic structure from the base of the convection zone to about 30 
Mm under the photosphere and that emergence higher up is not possible in these global anelastic 
simulations. Contact with three-diemnsional compressible local simulations of emerging flux in the 
solar atmosphere ( Archontis fc Hoodl|2010 ) should now be considered to have a full picture of the 
MHD evolution of large-scale magnetic regions. 

When convection is considered, the emerging radial field may lose the typical characteristics 
of a well-defined bipolar structure that we get in the isentropic cases. Indeed, in Case ITwP where 
the evolution was dominated by the downflows and upflows of the surrounding convection, the 
emerging structure is very complex and regions of mixed-polarity are self-consistently formed. The 
large variety of large-scale emerging regions at the solar photosphere may then well be the results 
of magnetic fi-loops significantly processed by convection, even if robust statistical properties such 
as tilt angle or asymmetry between polarities indicate that these loops are still able to keep their 
coherence during their rise. Another possibility leading to the emergence of complex active regions 
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in the Sun could be that such buoyant structures might interact in the bulk of the convection zone. 
Depending on which kind of reconnection could occur between the fi-loops, either complex region 



with mixed polarity and mixed helicity could emerge (such as the one observed by Chandra et 



al. 2010| ), or the magnetic field could completely change its trajectory, leading to no emergence 



at all (see Linton &; Antiochos 2005 for a discussion of reconnection between twisted flux tubes). 



This is the analysis we intend to focus on in the near future. Further no background dynamo field 



was present in the simulation. Pinto &; Brun (2012) have started to address the influence of such 



a dynamo field on a uniformly buoyant magnetic rope. We intend to carry on this analysis for 
fi-loops. 
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